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Abstract 

An efficient approach to tensor perturbation calculations by proper use of computer al- 
gebra methods is described, reaching the sufficient generality required for a comprehensive 
analysis of the Schwarzschild and Reissner-Nordstr0m metric stability problem. 



1 Introduction 



The stability problem for a Schwarzschild black hole in the form of a "pure metric" pertur- 
bation analysis was settled by T. Regge and J. A. Wheeler in a classic 1957 article |jl|. The 
main practical achievement of this work was undoubtedly the formulation of a gauge transfor- 
mation approach that allows a complete radial/angular separation of the Einstein equations 
in the two cases of odd and even parity, preliminarily established, reaching a formal solution 
in the axial (or magnetic) case: the so-called Regge- Wheeler Equation. Due to mathematical 
complications, however, the full analysis was only completed thirteen years later, after the 
work of Mathews |^], Edelstein-Vishveshwara and Zerilli Q that either provided a more 
rigorous approach to the use of tensor harmonics, or resolved some compatibility problems in 
the analytic treatment of the system, or else provided the final form of the equation for the 
radial perturbation functions in the polar (or electric) case: the so-called Zerilli Equation. 
The search for simplicity also led these authors to exploit some useful but not general rela- 
tions between curvature tensors, like those derived by Eisenhart |||], valid to the first order 
and/or only in the Schwarzschild case (latin tensor indices are used for consistency with the 
implemented algorithms) : 

SGmn = ^Rmn (1) 

^Rmn = ^r^^.p — (5r^p.^ (2) 
^^]k = ^Q^^'ihjp-k + hkp-j — hjk-p) (3) 

where h is the perturbation tensor and 6T, 5R, 6G are the perturbed parts of the affine 
connections, Ricci and Einstein tensors. 

On this way, after having performed the gauge transformations for each parity case, we are 
left with two systems of, respectively, three and six independent radial ordinary differential 



equations in two (namely hQ{r), hi{r)) and three (namely H{r), Hi{r), K{r)) arbitrary 
perturbation functions, to be determined. The magnetic system is straightforwardly reduced 
to a single first-order equation in /ii(r) which, by a simple variable substitution, leads to the 
final result: 

Q"{r*) + - V{r)]Q{r) = (4) 

(Regge- Wheeler Equation), where r* and Q (the first called the tortoise coordinate by Wheeler 
as a citation of the Zeno paradox) are (implicitly or not) defined by: 

d/dr* = a{r)d/dr Q{r) = b{r)hi{r) (5) 

where a(r), 6(r) are arbitrary functions to be determined in each particular case and V{r), 
in this Schrodinger-like equation, plays the role of an effective potential. 
Finally, the electric system, by a more complex change of variables procedure, necessary to 
deal with the terms in the wave number k, was found by Zerilli |^] to be represented by an 
equation formally equal to (Q) (the Zerilli Equation) only with a different (but still algebraic) 
expression of the potential. Since all the variables of the system are mutually expressed by 
regular algebraic relations, the stability problem, analyzed by substituting different forms 
of k into (Q), can be extended in its validity to the whole perturbation. The aim of the 
present work is to show that, with the essential help of computer algebra software, a similar 
analysis can be carried out for both parity cases, dealing with the more general spacetime of a 
spherical, non-rotating, eventually charged collapsed object, therefore allowing specialization 
not only to the Schwarzschild but also to the Reissner-Nordstr0m metric. 

2 The static spherically-symmetric system 

To begin a less rigid analysis than that induced by the formulation of equations (1-3), we 
deal first of all with the full expression of the Einstein tensor which, viewed as a function of 
the metric and its ordinary partial derivatives, reads as: 

ln[2] = = Sln^lify [Tsliiplif¥[3lff ineToMe trie [FicciToBf fine [EicciR [In, In]]] - 
1/ 2 Diim[]lf fineTaHjetricCScalarRtDAf finie[ScalarR]] ItetricgCliii, In]]]] 
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2 (gpq,rj - gpj ,qT ' ^qr^pj + Sr; ,pq) ) ) - 
2 Cg") {{<i") CCgp;,!,) (g,r,ivi) + 2 (gpj,,;) C- (gqm^r) + grm^q) - Cgr;,p) (g,rv,,ii+ gqn,m -gmii,q) ) + 
^ C r fi'^ri ~ 9^pm ,■ qn ~ Q^qn m ^rvm ^p q ) ) ) 

Here, to the unperturbed metric, which refers to the usual covariant expression of the line 
element belonging to a generic spherically-symmetric expression depending on the radial 
arbitrary functions A(r), u{r) : 

ds^ = e^^'Ur^ + r\de^ + sin^ 9 dcj)^) - e^'Ut" 
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must be added the perturbation terms, whose time dependency is represented only by a e~ 
factor while the angular one is limited to functions of the polar angle 9. 

Since " Components" , MathTensor's routine for calculating curvature tensors from an input 
metric, has the useful feature of applying arbitrary combinations of Mathematica commands 
to each component 0, the first order form of the Einstein tensor can be readily obtained by 
appending to the metric input file, where the Regge- Wheeler-gauge perturbation terms are 
added with a "small parameter" g as a factor, a line like: 

CompSimp [a_] : =Siniplif y [Normal [Series [Expand [a] , {q,0 , l}] ] / . q->l] 



2.1 The magnetic case 

There are only three non-zero components {(r, i;^), (6*, (j)), (</), t)} of the tensor equation that 
replaces the (1): 

SGmn = G'^nn — G^^^ = (6) 

(where the subtracted quantity is the unperturbed metric tensor) and that substantially 
reproduce, in their radial form, the solution system [^-(sys. 2), once the angular terms 
coming from the perturbation, where they are represented by f{0) = s'm9 dPL{cos 6)/d6 {Pl 
being the Legendre polynomial to the multipolar order L), are fully simplified through the 
following relations: 

f{9) = L [-Pl_i (cos e) + L cos 9 Pl (cos 9)] 
f'{9) = -L{L + 1) sin 9 Pl(cos 9) 

f"{9) = L{L + 1)[L Pl^i{cos9) - (L + 1) cos 9 Pl {cos 9)] 
and, being completely factored, are consequently eliminated. 

The first two equations (from now on, almost everywhere in the rest of this paper, the no- 
tation A = (L — 1)(L -|- 2)/2 will be adopted -not to be confused with the definition of the 
metric's radial function A(r)!) turn out to be of the first order with non-mixed dependence 
on the functions' derivatives in such a way that a single second-order differential equation can 
be obtained by: 

ln[12]:= Eql = 41 E"'''^"^ k r ho [r] + 

(4-4 U + 1> +2E"^^^ k^r^ +E-^^^^r ( (A'[r] - v'[r] ) (2 + r v'[r]) - 2 r v" [r] ) ) hi[r] - 

2IE"'''^''^kr^ li^[r] ==0; 
Eq2 = -2IE^'^''^ hlio[r] + E"'^''^ (A'[e] -v'[r]) hi[r] - 2 E'''^'^^ hi[r] == 0; 
s= Flatten[SolvE[Eq2 , lio[r]]]; 
ds = 0^ s ; 

p= Espaiul[First[FiaiSinpli£y[Eql / . Flatten[{s, ds}]]]]; 
c£ = Coe££icient[p, hj[r]]; c£ * 

E^ = Collect[E»paiid[Sln¥li£Y[p/c£]], {hi[r], hi[E] , hi[r]}, FiQlSlnplifY] ==0 

Out[i:]= 2 E"^'"'' E% 



/ 2 3X'[z] 3 V' h: \ , , r , 1 r , 

Out[i3]= + hi[E]+ -Chi[E] 

^ r 2 2 ' 2-£} 

C2E-'[''' CE-'^I''' - 21) + E (eA'[e]^ - 4v'[E] +A'[e] (4- e v'[e]) - e + V[e])))) + 



h£[E] 
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Analyzing the dependence on k, it is readily seen that to impose the equation (^) with the 
conditions (||), a procedure of polynomial coefficient identification is required, carried out in 
the following way: 

ln[14]- Eq = E:^ai«i[e,((cl,(b[r]hi[r]na[r])a[r]] - V[r ] ) b [r] hi[r] ; 

SirijEq = Collect [Eg/ Coeff icient [Eg, hJEr]], {hiEr] , hi[r], h^[r]}, Togetlier] == 

Cb[E] a'[E] + 2a[E] b'[E]) hi[El 



Out [14]= 



a[E] 13[E] 

-IVFrl ^-Fi\r^ flTrlViTrl -i-arT-l^ViTrn 

+ h^[E] 



hi[E] Ck^l3[El - l3[E]V[E] + a[E] a'[E]b'[E] +a[E]^l3"[E]) 



a[E]!l3[E] 



ln[lS]- Tort = Solve[Coef ficient [Coefficient [Fir St [SiriiEg] , hi[r]], k^] == 
Coefficient[Coef ficient [First [E^] , hi[r]], k^], 
a[r]] 

Ou.[1i]= {{a[E]^-ET^-^['l-''['l^}, {a[E] ^ET ^-^t'l+^t'l^}} 

ln[l6]:= a[r_] = Last [Tort [[ 2 , 1]]]; CoeffHl = 

DSolvB [Coefficient [First [SlniiEq] , Iii[r]] == Coefficient [First [E^] , hi[r]], l)[r], r] 

l-l- VEi:[l]-A[r]+y[r] JEC[l]-;L[rl+V[r] 
Out[16]= ||b[E] |, |b[E] ^ II 

ln[17]:= b[r_] = Last [CoeffHl [[2, 1]]] /. C[l] -s- 0; FiaiSinplifY[ 

Solve [Coefficient [First [E^] , hi[r]] Coefficient [First [Sii*Eg] , hi[r]], V[r]]] 

E I "^1 (4 + 4E-^["^1 X + r (- (J.' [e] - V [e] ) (1 + E v'[e]) + 2 e v"[e])) 
Out[17]= ||V[E] ^ — II 



2.2 The electric case 

Here we deal with the seven non-zero components of the perturbed Einstein tensor (the 
complementary set with respect to the previous three "magnetic" components), whose angular 
parts, coming from the perturbing gauge function F{9) = Pl{cos9), are factored, once the 
substitutions 

F'iO) = [-PL-iicosO) + cose PLicosO)] 

F"{e) = PL-i{cose) - L (cot2 ^ + L + 1) Pl{cos9) 

sin 6 

are performed, the resulting expressions simplified and the original four radial perturbation 
functions are reduced to three by the identification HQ{r) = H2{r) = H{r). 
Of the six independent linear differential equations so obtained (the diagonal 9 and (f) terms 
being equal), four turn out to be of the first order, unlike the Schwarzschild-specialized system 
which instead has three first order equations, and a completely algebraic variable-elimination 
procedure allows the derivation of a first integral condition which is the generalization of 
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|^]-(eq. 10). As in Zerilli's procedure more elaborate but always tied to a polynomial co- 
efficient identification principle, this algebraic condition plus three of the previous equations 
can be treated, with a double function substitution and the change of the radial coordinate, 
to form a new system of seven differential equations in four variables (plus the derivative of 
the new radial coordinate with respect to r), this time non-linear but independent of /c, to 
which the formal definition of the effective potential must be added. Three of these unknown 
functions are then found, after a cascade of algebraic eliminations (which don't show the 
residual arbitrariness found by Zerilli in the Schwarzschild treatment), to be quite simply de- 
pendent on the fourth, which satisfies a final second-order very complex equation, fortunately 
analytically solvable when specialized to the Schwarzschild and Reissner-Nordstr0m metrics. 

3 Results of the analysis 

A comprehensive table of functions and variables in the three cases - general static spherically- 
symmetric form (G), Schwarzschild (S) and Reissner-Nordstr0m (RN)- can be sketched, the 
latter two obtained, respectively, by the two substitutions {A(r) = — ln(l — 2m/r) ; v{r) = 
-A(r)} and {A(r) = - In (1 - 2m/r + Q^/r^) ; y{r) = -A(r)} : 



TT 



Schrodinger wave-like function Tortoise coordinate 



G 



mag. 



Q(r) = exp [\{u{r) - A(r))] ^ r* = / exp [\{X{r) - u{r))] dr 



el. K{K{r), 



Hi(r)- 



A(r) 

e 2 



2Ae^('')-rV(r)+2n/'(r) 



dr 



S mag. Q(r) = (l-^)^ 



r* = r + 2m ln(r — 2m) 



el. K{K{r), 



r 



RN mag. Q(r) = 1 - ^ + 



2m _^ \ fei(r) 



el. K{K{r), 



gi(r-) - 
r 



r* = r + mln(r'^ — 2mr + Q^) + r]{r) 



rj{r) 



2m^ 



: arctan 



and the correspondent couples of specialized expressions of the Regge- Wheeler and Zerilli's 
potentials are: 



Schwarzschild: 



>, , , 2m \ / A + 1 3m 



Viir) = 2 



^ 2m\ A2(A + l)r^ + SA^mr^ + OAm^r + 9m^ 
r J (Ar + 3m)^r^ 
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2m Q2\ p(^) 



r r2 J 4 [r(Ar + 3m) - 3Q2]2 ^4 



with 



P(r) = 8A2(A + l)r6 + 2A\^mr^ + 2A [SGm^ - 7(2A - 3)Q2] ^4 
+ 12m [6m^ - (19A - 3)Q2] - SQ^ [108m^ - (38A - ?>)Q^\ 
+ 342mQ^r - 117Q^ 



As a final remark, it is straightforward to verify that the two different kinds of tortoise 
coordinates in the general case reduce to one for all metrics having A(r) = —v{r) and that as 
expected, limg^o V^^ {r) = V^{r) holds for both the magnetic and the electric parities. 
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